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ABSTRACT 

We present the first HST long-slit spectrum of a gaseous disk around a candidate 
super-massive black-hole. The results of this study on the kinematics of the gaseous 
disk in M87 are a considerable improvement in both spatial resolution and accuracy 
over previous observations and requires a projected mass of MBH(sinz)^ = (2.0 ± 
O.5)xlO®M0 (Mbh = 3.2xlO®M0 for a disk inclination i = 52°) concentrated within 
a sphere whose radius is less than O'.'OS (3.5 pc) to explain the observed rotation curve. 
The kinematics of the ionized gas is well described by a thin disk in keplerian motion. 
A lower limit to the mass-to-light ratio of this region is M/Ly —110, significantly 
strengthening the claim that this mass is due to the presence of a central black-hole 
in M87. 
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1 INTRODUCTION 

One of the cornerstones of the contemporary view of the 
physics of active galactic nuclei is that their energy out¬ 
put is generated by accretion of material on to a massive 
black-hole (e.g. Blandford 1991, Antonucci 1993 and ref¬ 
erences therein). Considerable controversy has surrounded 
attempts to verify the existence of black-holes in nearby gi¬ 
ant elliptical galaxies using ground-based stellar dynamical 
studies (Sargent et al. 1978, Young et al. 1978, Duncan & 
Wheeler 1980, Binney and Mamon 1982). To-date the best 
available data remains inconclusive, largely because of the 
difficulty of detecting the high-velocity wings on the absorp¬ 
tion lines, which are the hallmark of a black-hole (van der 
Marel 1994, Dressier & Richstone 1990, Jarvis & Pelletier 
1991, Kormendy & Richstone 1995). One of the major goals 
of HST has been to establish or refute the existence of black- 
holes in active galaxies by probing the dynamics of AGN at 
much smaller radii than can be achieved from the ground. 
Already persuasive evidence has been found that apparently 
quiescent nearby galaxies contain black-holes (van der Marel 
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1997a, Kormendy et al. 1996) but, locally, there are few can¬ 
didate galaxies which might harbour a super-massive black- 
hole large enough to have sustained Quasar type activity 
earlier in their evolution (Rees 1997). M87 is the nearest and 
brightest elliptical believed to harbour such a super-massive 
black-hole. 


M87, the dominant giant elliptical galaxy at the center 
of the Virgo cluster, was the first galaxy for which tenable 
stellar dynamical and photometric evidence was advanced 
for the presence of a supermassive black-hole (Sargent et al. 
1978, Young et al. 1978). Subsequent stellar dynamical mod¬ 
els showed however that the projected stellar density and the 
observed rise in velocity dispersion did not necessarily imply 
the presence of a black-hole, but could still be explained by 
the presence of a stellar core with an anisotropic velocity 
tensor(Duncan & Wheeler 1980, Young 1980, Binney and 
Mamon 1982). 


An important step in the quest to verify the existence 
of a massive black-hole in M87 has been the discovery from 
HST emissio n line (Crane et al. 1993, Ford et al. 1994) and 
continuum (Macchetto 1996) imagery of a circum-nuclear 
disk of ionized gas which is orie nted approximately perpen¬ 
dicularly to the synchrotron jet (Sparks et al. 1996). Similar 
gaseous disks have also been found in the nuclei of a num¬ 
ber of other massive galaxies (Jaffe et al. 1993, Ferrarese et 
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al. 1996). Because of surface brightness limitations on stel¬ 
lar dynamical studies at HST resolutions, the kinematics of 
such disks are in practice likely to be the only way to deter¬ 
mine if a central black-hole exists in all but the very nearest 
galaxies. 

Previous HST spectroscopic observations (Harms et al. 
1994, Ford et al. 1996) of the gas disk at several discrete 
locations on opposite sides of the nucleus showed a velocity 
difference of >1000 km s“^. On the assumption that these 
motions arise in a thin rotating keplerian disk, this lead to 
a central dark mass of M87 in the range 1 to 3.5x10®Mq. 
An important shortcoming of this estimate is the assump¬ 
tion about the inclination of the disk with respect to the line 
of sight which cannot be properly determined neither by us¬ 
ing the few velocities derived from the FOS observations nor 
from the WFPC2 imaging data. Harms et al. (1994) assumed 
a disk inclination of 42° ± 5°, as derived by an ellipse htting 
of the Ha-|-[NH] image, and this resulted in a misleadingly 
accurate mass value. Actually, such value for the inclination 
was determined at distances between C/.'S and Cf.'S from the 
nucleus, because in the inner regions the disk emission was 
subsumed by the bright central point-sonrce. Indeed Ford et 
al. (1996) quoted a larger uncertainty in their mass determi¬ 
nation but it is not clear from that short paper whether they 
have fully taken into account all of the possible values of the 
disk inclination which are consistent with the data. Implicit 
in this measurement of the mass of the central object is 
also the assumption that the gas motions in the innermost 
regions reflect keplerian rotation and not the effects of non- 
gravitational forces such as interactions with the jet which 
are known to dominate the gas motions in the inner regions 
of many AGN (Whittle et al. 1988, Axon et al. 1989). Es¬ 
tablishing the detailed kinematics of the disk is therefore 
critical. Currently their result still provides the most con¬ 
vincing observational evidence in favour of the black-hole 
model. 

In order to determine whether the motions observed in 
M87 are due to keplerian rotation, we have re-investigated 
the velocity field at high spatial resolution with HST. Our 
results demonstrate that both the observed rotation curve 
and line proHles are consistent with a thin-disk in keplerian 
motion and allow us to determine an improved estimate for 
the mass of the black-hole. On the assumption that the dis¬ 
tance to M87 is 15 Mpc, 0”l~7pc. 


2 OBSERVATIONS AND DATA REDUCTION 

The circumnuclear disk of M87 was observed on 1996 July 
25*^, with the COSTAR corrected Faint Object Camera f/48 
long-slit spectrograph on board the Hubble Space Telescope. 
An F305LP filter, was used to isolate the first order spec¬ 
trum which covers the 3650-5470 A region. All the spectra 
have a pixel size of 0'.'0287 in the spatial direction and cover 
the spectral lines of [O H]AA3726,3729 A, H/3A4861 A, and 
[O HI] A4959,5007 A at a resolution of 1.78 A per pixel and 
were taken in the 1024x512 non-zoomed format of the Faint 
Object Camera. An interactive acquisition image (F140W 
Hlter) was first obtained to accurately locate the nucleus. 
The slit, with dimensions of 13”5xCf.'063 at a position angle 
of 47°, was then stepped across the nucleus at three positions 
separated by 0”2. The integration times for the inner spec- 



Figure 1. The slit positions discussed in the text are marked 
with solid lines on an Ha-|-[NI1] image from the HST archive. 
North is up and east is left. 


trum was 7761 second while those of the two outer spectra 
where 2261 seconds. Comparing the continuum flux observed 
in the spectra with the luminosity distribution derived from 
a F342W FOC, f/96 archival image we could estimate that 
the inner slit position was located to within 0”07±iy.'01 of 
the nucleus. The slit positions used are shown in Figure ^ 
overlaid on an Ha-|-[NH] WFPC2 image of the gaseous disk 
of M87 obtained from the HST data archive. 

The distortion induced by the optics and the magnetic 
focusing of the detector was corrected by using the equally 
spaced grid of r eseau marks, etched onto the pho tocathode 
of the detector (Nota, Jedrzejewski & Hack, 1995). The dis¬ 
tortion induced by the spectroscopic mirror and the grating 
were corrected using spectra of the planetary nebula NGC 
6543 and of the globnlar cluster 47 Tucanae for the slit and 
dispersion directions respectively which are then character¬ 
ized by maximum uncertainties of 0.45 pix and 0.28 pix (i.e. 
0.8 A and 8 mas) when comparing measurements all across 
the detector. After the background subtraction, the only line 
with enough signal to be suitable for velocity measurements 
was [O H] A3726,3729 A, which was then fitted (row by row) 
using the task LONGSLIT in t he TWODSPEC FIGARO 
package (Wilkins & Axon 1992). In all cases the line pro- 
Hle is well represented by a gaussian function. Full details of 
the observation st rategy and data reduc tion are given in a 
companion paper (Macchetto et al. 1997). 


3 RESULTS 

The measured variations of radial velocity and continuum 
flux for the central slit are shown in Figure ^ It is immedi¬ 
ately apparent that the observed velocity field does not show 
the very steep rise in velocity one would naively expect to 
be the signature of a massive black-hole. Instead over the 
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Figure 2. Observed rotation curve from the [O II]AA3726,3729 
A doublet for the spectrum with the smallest impact parameter 
to the nucleus of M87 compared with the continuum flux distri¬ 
bution along the slit. The uncertainties in the measured velocities 
and positions are 65 km s“^, at 3727A, and 0.28 pixels (8 mas) 
respectively. 


central ±0”2 the velocity varies approximately linearly with 
position, with an amplitude of ±600 km s“^, while at larger 
radii it flattens, and eventually starts to turn over. As shown 
in the following section, this apparent flattening is a conse¬ 
quence of both the finite distance between the slit and the 
nucleus and the smearing of the instrumental spatial PSF. 

The comparison between the new FOC data and the 
archival FOS observations is shown in Fig. ^ where we in¬ 
clude only those FOS apertures which overlapped our slit 
at NUC. Within their substantially larger uncertainties, the 
previous velocity measurements (Harms et al. 1994, Ford et 
al. 1996) are in reasonable accord with our results. The com¬ 
parison shows clearly the increase in spatial resolution and 
sampling of the FOC observations with respect to the FOS 
ones. Moreover, our data provide considerably more reliable 
relative spatial positions, since all the data points were ob¬ 
tained simultaneously on the detector. On the contrary, the 
limitations of the FOS targ et acquisition procedure can pro¬ 
duce errors as large as 0”l (van der Marel et al. 1997b) and 
such uncertainty in the position is crucial in the inner 0”2 
when comparing observations taken at different epochs. 


4 DISCUSSION 

To understand the implications of our results for the mass 
distribution in the center of M87 we built models of the 
gas kinematics both under the assumption of the existence 
of a central black-hole, and an extended mass distribution 
following a Plummer Potential. 

In the hypothesis of a thin disk rotating in circular or¬ 
bits, neglecting the slit width, instrumental PSF and lumi- 

© 1997 RAS, MNRAS 000, §-| 


nosity distribution of the line, the observed velocity is given 
by: 

V = Vsys-HRf-^^smi (1) 

H 

where 


X = —b sin 0 ± s cos 6 Y = b cos 9 -\- s sin 9 


and 


7?=(x^ + y2)0.5 

b is the distance of the slit from the nucleus, s is the coordi¬ 
nate along the slit and R represents the “true” distance of 
each point on the slit from the nucleus (s = 0 when R — b). 
Vsys is the systemic velocity, i is the inclination of the disk 
with respect to the line of sight {i = 90° for the edge-on 
case), 9 is the angle between the slit and the line of nodes. 
<E>(i?) is the gravitational potential and results 


HR) 

HR) 


GMbh 


Keplerian 


GM 

(i?2 ±i?^)° ® 


Plummer 


where Mbh and Rc are the black-hole mass and core radius, 
respectively. 

The effects of the f/48 spatial PSF and the finite slit 
size are taken into account by averaging the velocity field 
using the luminosity distribution and PSF as weights. The 
model rotation curve is thus given by: 


... fs-^s dbf f^: db'ds'V{s',b')I{s',b')P 

Vps[o)~ ^SHAS . nB^h .. n nHoc, .... 


/r b')p 


where V{s',b') is the keplerian velocity derived in eq. 
I{s',b') is the intrinsic luminosity distribution of the line, 
P = P{s' — s,b' — b) is the spatial PSF of the f/48 relay along 
the slit direction. B is the impact parameter (measured at 
the center of the slit) and 2h is the slit size, S is the position 
along the slit at which the velocity is computed and 2AS' is 
the pixel size of the f/48 relay. For the PSF we have assumed 
a gaussian with 0'.'08 FWHM i.e. 


P = 



exp 



Hb'-b)^ \ 

2 J 


(3) 


Because the intrinsic surface brightness distribution of 
the emission line disk interior to the HST PSF is not known, 
we modeled this with two extreme parameterized forms, ei¬ 
ther as a power law or as an exponential. We then fitted the 
model rotation curves using a minimization. The free pa¬ 
rameters of the fit, in the case of the black-hole model, were 
the inclination of the disk with respect to the line of sight, 
the angle between the disk and the line of nodes (positive 
angles indicate that the line of nodes has a larger position 
angle than the slit), the position along the slit of the point 
closest to the nucleus, the impact parameter of the slit, the 
systemic velocity and the mass of the black-hole. The uncer¬ 
tainties are dominated by the poorly determined inclination 
of the gas disk within the HST PSF. Taking into account 
the possible ranges of variation of the model free parame¬ 
ters, which were derived using a large Montecarlo generated 
grid of solutions, we find a best estimate of the projected 
mass of the black-hole of Mbh (sin f)^ = (2.0 ± 0.5) x lO^M©. 
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Figure 3. Comparison between the FOC rotation curve and the 
archival FOS observations. The labels “pos#” refer to the nota¬ 
tion used by Ford et al. (1996). “a” and “b” indicate positions 
which were nominally coincident but resulted slightly shifted due 
to acquisition uncertainties. The error-bars in position represent 
the aperture diameters, those in velocity represent the disper¬ 
sion of the measurements on different lines. The solid, dotted and 
dashed lines are model fits which are discussed in the text. 


Three representative model fits are shown in Fig. The 
solid, dashed and dotted lines represent three black-hole 
model fits which are compatible with the observations and 
differ in the value of the angle between the slit and the line 
of nodes of the gaseous disk (9 ~ —9°, 7° and 1°, respec¬ 
tively). The best fit model, represented by the solid line in 
the upper panel has a disk inclination, i = 51° (39°; 65°), an¬ 
gle of line of nodes, 6 = —9°(—15°; 13°), impact parameter, 
b = (y.'08(0”06; C/.'OS) and a systemic velocity, Vsys = 1290 
km s“^ (1080; 1355), leading to a mass estimate for the 
black-hole of Mbh = (3.2 ± 0.9) xIO^M q. The values in 
parentheses represent the allowed errors on the values of 
the parameters. As shown by this analysis, the main source 
of uncertainty on the deprojected mass is the value of the 
disk inclination which more than doubles the total error. 

In the case of an extended mass (Plummer) distribution 
the observations are reproduced only if at least 60% of the 
total mass is confined within a sphere of 0”05 (~3.5pc). 
Since the maximum spatial extent is determined by the PSF 
FWHM (0V08) this possible mass distribution is therefore 
indistinguishable from a point mass. 

We have demonstrated that the data are fully compat¬ 
ible with a thin disk in keplerian rotation around a central 
mass condensation, at least within the inner 0"2-0”3. At 
larger distances, the signal-to-noise in our data decreases 
and, consequently, the error bars increase. The few points 
that appear to deviate from the keplerian rotation curve, do 
not provide enough information or weight to warrant the fit 
of more complex models, such as thick disks or sub-keplerian 
disks with outflows. 


Can we rule out the alternative explanation that the 
observed rotation is due to an extended mass? A strong 
argument against this interpretation is that a mass of 
(3.2 ± O.9)xlO®M0 distributed within 3.5pc yields a mean 
density of ~ 2x10^Mq pc“®, higher than the highest den¬ 
sity encountered in the collapsed cores of galactic globular 
clusters. Moreover, the V band flux in a region 16x 16 pc^ of 
the nucleus yields a mass-to-light ratio M/Lv~110 Mq /Lvq 
(Lv0=O.113 Lq) which is far higher tha n those predict ed for 
stellar clusters by evolution synthesis ( Bruzual 1995 ). The 
mass condensation in the nucleus of M87 can therefore not 
be a cluster of normally evolved stars. Even more exotic pos¬ 
sibilities (e.g. massive clusters of neutron stars, other dark 
objects etc.) have been discussed by van der Marel et al. 
(1997a) in the case of M32 and found to be equally improb¬ 
able. Therefore the most natural explanation for the central 
mass condensation of M87 is that it is due to the presence 
of a black-hole. 

The data and the analysis presented by Harms et al. 
(1994) and Ford et al. (1996) indicated the presence of a 
velocity shear which could be consistent which the presence 
of a black-hole. Here we have demonstrated the existence 
of a keplerian rotation curve in the disk and we have also 
excluded the possibility that it might be due to a mass con¬ 
centration more extended than our spatial resolution. 


5 CONCLUSIONS 

We have presented the first HST long-slit spectrum of a 
gaseous disk around a candidate supermassive black-hole. 
We have obtained a rotation curve in the [O H]A3736,3729 
which extends up to ~l” from the nucleus. We have mod¬ 
eled the rotation curve in the case of a thin disk in circular 
orbits around a supermassive black-hole and we have de¬ 
rived a projected mass of Mbh (sin i)^ = (2.0 ± 0.5)x IO^Mq 
(Mbh = 3.2xlO®M0 for a disk inclination, i = 52°) con¬ 
centrated within a sphere whose radius is less than 0?05 
(3.5 pc). We have shown that the observed characteristic of 
the data are well explained under the working hypothesis 
of circular keplerian orbits and that there is no substantial 
contribution from a mass distribution more extended than 
our spatial resolution. Given the inferred mass-to-light ra¬ 
tio M/Lv~110 M 0 /LV© the most natural and likely expla¬ 
nation is that of a supermassive black-hole. In conclusion, 
with respect to previous determinations, we have improved 
the accuracy of the mass estimate, demonstrated its reliabil¬ 
ity by verifying the assumption of the thin disk in circular, 
keplerian rotation and excluded the possibility of a mass 
distribution more extended than the spatial PSF. 

To make further progress there are a number of possi¬ 
bilities the easiest of which is to make a more comprehensive 
and higher signal-to-noise 2D velocity map of the disk to bet¬ 
ter constrain its parameters. The biggest limitation of the 
present data is that, even by observing with HST at close to 
its optimal resolution at visible wavelengths, some of the im¬ 
portant features of the disk kinematics are subsumed by the 
central PSF. Until a larger space based telescope becomes 
available, the best we can do is to study the gas disk in Lya 
and gain the Rayleigh advantage in resolution by moving 
to the UV. This may be the only way to proceed because 
of the difhculty of detecting the high velocity wings which 
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characterize the stellar absorption lines in the presence of a 
supermassive black-hole. 
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